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ABSTRACT 

We present a new Bayesian approach to constrain the intrinsic parameters (stehar mass, age) of 
the ecUpsing binary system ~CEP0227- in the Large Magellanic Cloud. We computed several sets of 
evolutionary models covering a broad range in chemical compositions and in stellar mass. Indepen- 
dent sets of models were also constructed either by neglecting or by including a moderate convective 
core overshooting (/3ot,=0.2) during central hydrogen burning phases. Sets of models were also con- 
structed either by neglecting or by assuming a canonical (?7=0.4,0.8) or an enhanced (?7=4) mass loss 
rate. The most probable solutions were computed in three different planes: luminosity-temperature, 
mass-radius and gravity-temperature. By using the Bayes Factor, we found that the most probable 
solutions were obtained in the gravity-temperature plane with a Gaussian mass prior distribution. 
The evolutionary models constructed by assuming a moderate convective core overshooting (/3ou=0.2) 
and a canonical mass loss rate (7^=0.4) give stellar masses for the primary (Cepheid) -M=4.14lg Q5 
Mq- and for the secondary -M=4.15to o5 ■'^0" ^^8.t agree at the 1% level with dynamical measure- 
ments. Moreover, we found ages for the two components and for the combined system -t=151^3 Myr- 
that agree at the 5% level. The solutions based on evolutionary models that neglect the mass loss 
attain similar parameters, while those ones based on models that either account for an enhanced mass 
loss or neglect convective core overshooting have lower Bayes Factors and larger confidence intervals. 
The dependence on the mass loss rate might be the consequence of the crude approximation we use 
to mimic this phenomenon. 

By using the isochrone of the most probable solution and a Gaussian prior on the LMC distance, we 
found a true distance modulus -18.53tQ;g2 mag~ and a reddening value -E(B-V)= 0.142j;g Q5Q mag- 
that agree quite well with similar estimates in the literature. 

Subject headings: stars: classical Cepheids — stars: evolution — stars: oscillations — stars: funda- 
mental parameters 



1. INTRODUCTION 

Classical Cepheids are very popular in the recent lit- 
erature, since they play a crucial role to address several 
open astrophysical problems. The reason is threefold. 
i)- They are the most popular primary distance indica- 
tors, since they are bright and can be easily identified. 
a)- They are excellent tracers of intermediate- mass stars, 
since their evolutionary status is well defined (central he- 
lium burning, blue loops), in)- They are fundamental 
laboratories to constrain the micro (equation of state, 
opacity, cross sections) and macro (mass loss, mixing, 
convective transport) physics adopted in both evolution- 
ary and pulsation models. 
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Dating back to the seminal inves tigation by 
iHofmei ster. Kippenhahn fc Wcigcrt' (1964', and refer- 
ences therein) and by Becker & Ibcn (1979), the evolu- 
tionary properties of intermediat e-mass have been inves- 
tigated in a series of theoretical (IStothers fc ChirJ 119781: 
Chiosi fc Macdcr 1986; Castellani. Chieffi fc Stranierc 
1990; Lan genil991c iBono et al. 2000; Macder fc MevneJ 
2000; Me vnet fc Maedei) 12002: Karakas fc Lattanzic 
20031: JPietrinfcrni et al. 2004, 2006 ; IValle et aT 



2008; 



2011 



Charbonncl fc Lagardc 20X0; Cassisi fc Salaris 



Brott et al 



_ ■ .Neilson. Cantiello. fc LaiigeH 120111. 
2011allH) and empirica l (lAlcock et al.l 119991: iTesta et al 
19991: nBrocato et all 120031: iHunter et al.l | 2009all 



Bonanos e t all 120091 120101: ISanna et al 



200 



Efrcmova et al.ll2011f ) investigations 

One of the most interesting problem concerning evo- 
lutionary and pulsation properties of classical Cepheids 
is the so-called mass discre pancy problein. Du ring the 
late sixt ies it was noticed bv lChristvl ()1966l , fl970l) and by 
Stobi(i (|1969f ) that the evolutionary masses -estimated 
using the comparison between isochrones and observa- 
tions in the Color-Magnitude Diagram (CMD)- were al- 
most a factor of two larger than the pulsation masses 
-estimated using the Period-Mass-Radius relation- o f 
Galactic Cepheids (|Fricke. Stobie. fc Strittmatteiiri97lD . 
This conundrum was partially solved by 
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iMoskalik. Buchler. fc MaromI ()1992[ ) using the 
sets of radiative o paciti es released by the OPAL 
(jlglesias fc Roger j ll99lD and by the Opac- 
ity Project (jSeaton et alJ Il994| ) groups. How- 



ever, several investigations focussed on Galactic 



Bono et all 120011: fC aputo et al. 2005)[) and Magellanic 
.Beaulieu. Buchler. & Kollath 2001; Bono c t all 120021: 
IKeller fc WoodI 12002. ,2006i .Keller, .2006) Cepheids 
suggested that such a discrepancy was still of the order 
of 10-20%. A similar discrepancy was also found by 
lEvans et al.l (jTOOSa) using dynami cal mass estimates o f 
Galactic binary Cepheids and by iBrocato et al.l ()2003f ) 
using cluster Cepheids located in NGC 1866, a young 
Large Magellanic Cloud (LMC) cluster. 

The Cepheid mass discrepancy problem can be ad- 
dressed following two different paths. 

a) -He-core mass- current evolutionary predictions 
underestimate the He-core mass for intermediate-mass 
structures, and in turn their Mass-Luminosity (ML) ra- 
tio. This possible drawback has a substantial impact 
both on evolutionary and pulsation predictions, since the 
latter models assume an ML relation. 

b) -Envelope mass- Current Cepheid masses are 
smaller than their main sequence (MS) progenitors, 
because they have lost a fraction of their initial en- 
velope mass. The latter working hypothesis implies 
that the Cepheid mass discrepancy is intrinsic, i.e., it 
is not caused by limits in the physical assumptions 
adopted in constructing evolutionary and pulsation mod- 
els. The input physics adopted in constructing evolu- 
tionary models of intermediate-mass stars has already 
been addressed in sev eral papers (Chiosi & Macdcr 198(: 
Brocato fc Castellani 1993; Stothcrs fc Chin 1993, 199f 
Cassisi ct al. 2004; BcrtcUi ot al. 2009; Vallc ct al. 200$ 

Here, we briefly mention the most relevant ones affect- 
ing the He-core mass and the envelope mass. 

a) He- core mass 

i) -Extra-mixing- Several detailed studies based on 
the comparison between predicted CMDs and Luminos- 
ity Functions of NGC1866, reached opposite conclusions 
in favor (Ba-rmina. Girardi. fc Chiosi 20 021 and against 
()Testa et al.l Il999t IBrocato et al.l I2003D the occurrence 
of mild convective core overshooting. The need for a 
mild o vershooting was also suggested by IKeller fc WoodI 
(|2006') who investigated several young clusters in the 
LMC and in the Small Magell anic Cloud. M o re rec ently, 
it has also been suggested bv iCordier et al.l (|2002| ) that 
the degree of overshooting might also depend on the 
metal abundance, namely it increases when metal abun- 
dance decreases, ii) -Rotation- Evolutionary models 
constructed by accounting for the effects of rotation 
predict both an increase in Helium core mass and an 
enhancement in the surface abundanc e of both Helium 
and Nitrogen (jMevnet fc Maede"i)l2002D . It has also been 
suggested that rotation might account for the significant 
changes in surface chemical composi tions observed in 
Galactic and Magellanic supergiants (|Korn etTall 12005) : 
iHunter et al.l [2009ajb). Moreover, recent theoretical 
(|Maeder fc Mevnetl EoOOf ) and empirical (|Vennl Il999| l 
investigations indicate that the efhciency of such a mech- 
anism might depend on the initial metal abundance, 
i.e., it increases when metal abundance decreases, in) 
-Radiative Opacity- A new set of radiative opacities 



has been recently c omputed by the Opacity Project 
(|Badnell et al.ll2005l ). The difference between old and 
new opacities is at most of the order of 5-10% across 
the Z-bump (T« 250,000 K). To account for the mass 
discrepancy the increase in the opacity should be 
almost a factor of two. This indicates that the adopted 
radiative opacities have a marginal impact on the mass 
discrepancy problem. 

h) Envelope mass 
i) -Canonical Mass Loss- The total mass of an ac- 
tual Cepheid is also affected by a decrease in the en- 
velope mass. Evolutionary models accounting for the 
mass loss, during Hydrogen and Helium burnin g phases, 
by means of several semi-em pirical relations ()Reimera 
ll975tlNTeuwenhuiizen fc deJ agcr 1990) do not solve the 
Cepheid mass discrepancy problem. Plausible values 
for the free parameter -rj- give mass loss rates that 
are too small. This is not surprising, since current 
semi-empirical relations are only based on scaling ar- 
gume nts and they are not roo ted on a robust physical 
basis (Schroede r fc Cuntdl2005l ). Empirical mass loss es- 
timates based on NIR and ul traviolet em issio n for Galac- 
tic (IMcAlarv fc Wel ch T98g; lDeasvl [T988; Mar engo et al 



2010a) and Large Mage llanic Cloud (.Neilson fc Lester 
2008t Neilson et a"Lll2009D Cepheids cover a broad range 
(10"i° - 10""^ M© yr-i). ii) -Enhanced Mass Loss- 
More recent mid-infrared observations for 29 Galactic 
Cepheids collected with Spitzer indicate t he presence of 
exten ded emission around seven of them (|Barmbv et al.l 
1201 Ih . This finding together with the interferometric 
detection of circumstellar shells around several Galactic 
Cepheids fM crand ct al. 2006, 2007; KcrvcUa ct al. 200i, 
120081: iKervclla . Merand fc Gallennel l2009) further sup- 
ports the evidence that there is a wind associated with 
Cepheids possibly triggered by pulsation-driven mass 
loss. 

A new spin concerning the occurrence of circumstellar 
shells around Cepheids was recently given by the dctec- 
tion o f a large nebula around S Cephei (Marengo et al] 
l2010bD . the protype of Cepheid variables. The IR emis- 
sion was detected in several MIR bands and shows a ra- 
dial extent larger than 10** AU (5 arcmin). The nebula 
shows a parabolic shape and it is aligned with the di- 
rection of the motion, thus suggesting the possible oc- 
currence of a bow shock caused by the interaction be- 
tween the mass-losing Cepheid with the interstellar ma- 
terial. Current estimates indicate that S Cephei may be 
losing mass with a rate of ~5xl0~^ to ~6xl0~* Mq 
yr~^. This findi n g has been independently confirmed by 



iMatthews et al.l (|2011l ) using data collected in the ra- 
dio with the Extended Very Large Array. They found 
an outflow velocity of ^ 36 km/s and a mass loss rate 
of - lO-'^-lO"^ Mq yr-i. The modest dust content of 
the outflow indicates that the wind from 6 Cephei might 
be pulsation-driven instead of dust-driven as typical for 
other classes of evolved stars. 

However, the stepping stone concerning the evolution- 
ary and pulsation properties of classical Cepheids was the 
detection of a Cepheid in a well detached, double-lined 
eclipsing binary system in the Large Magellanic Cloud. 
The geometry and the precision of both photometric and 
spectroscopic data gave the opportunity to measure the 
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mass of the Cep heid (CEP0227) with th e unprecedented 
precision of 1% (jPietrzynski et all 120101 hereinafter Pa- 
per I). A similar precision was also attained for the clas- 
sical Cepheid in the LM C eclipsing binary CEP1812 by 
(jPietrzynski et aLll201lL hereinafter Paper II). 

This finding paved the road for new constraints 
on pulsation an d evol utionary masses. Soon after, 
ICassisi fc Salarii (|20T1 by using evolutionary models 
that account for a moderate convective core overshooting 
found evolutionary masses agree quite well with the dy- 
namical mass of CEP0227 both in the Radius- Age plane 
and in the Radius-Effective Tempera t ure p lane. More re- 
cently, iNeilson. Cantiello. fc Langeii ()2011l ) found that a 
combination of both moderate convective core overshoot- 
ing and pulsation-driven mass loss is required to solve the 
Cepheid mass discrepancy. 

This is the third paper of a series focussed on the evo- 
lutionary and pulsation properties of Cepheids. The new 
series present two novel approaches when compared with 
similar investigations available in the literature: i)- the 
intrinsic parameters will be estimated using a Bayesian 
approach in dealing with the comparison between the- 
ory and observations; ii)- the empirical measurements 
are based on a homogeneous approach for both pho- 
tometry and spectroscopy. In this investigation we pro- 
vide new evolutionary constraints on the binary Cepheid 
CEP0227. 

The plan of the paper is the following. In §2 we intro- 
duce the input physics adopted to construct evolution- 
ary models. We also discuss the physical assumptions 
adopted to deal with mass loss, mixing and convective 
transport. In this section, we also outline the assump- 
tions adopted to compute stellar isochrones. The de- 
scription of the method adopted to estimate the intrinsic 
parameters of the binary components is given in §3. In 
this section we discuss the different most probable solu- 
tions in the L-T, in the M-R and in the g-T planes. §4 
deals with the new empirical constraints on both stellar 
mass and age we obtained in the different planes for the 
binary system. In §5 we address the comparison between 
the most probable stellar isochrone and observations in 
the Color-Magnitude diagram. The true distance modu- 
lus and the reddening for the combined system are also 
discussed. In §6 we summarize the results and outline 
future possible avenues of the project. 

2. THE THEORETICAL FRAMEWORK 

2.1. Input physics 

The evolutionary tracks were computed with 
an updated version of the FRANEC evolution- 
ary 



code jP egl'Inn ocenti et al.l 120081 iValle et al.l 
2009"; TogneUi. Prada Mor oni, fc Degl'Innocentil 120111 : 
DeirOmodarme et al.. .2011 



For these models, the 
2006 rele ase of the OPAL e quation of state (EOS) was 
adopted ((Rogers et al.lll996l ). together with the radiative 
opacity tables released in 2005 by the same Livermore 
group (jlglesias fc Rogerd [T^96P I for \ogT[K] > 4.5 
and extended at lower temperatures w ith the radiative 
opacity tables bv lFerguson et al.l (I2OO50. Th e condu ctive 
opac ities are from iShternin fc YakovlevI ()2006l see 
also iPotekhinI ()1999D ). Ah the opacity tables were 

® http:/ /opalopacity.llnl.gov/opal.html. 



calculated for the lAsplund et al. I (|2009l hereafter As09) 
heavy-element solar mixture. The nuclear reaction rate s 
are from the NACRE compilation (jAngulo et al.lll99"9l) . 
except for t he ^^N(p,7)-'^^0 , which is from the LUNA col- 
laboration ( Imbriani et al.l [20051 : iBemmerer et al.l 120061 : 
iLemut et al1l2006( ) and '^'^C{a,jy^O fr omlHammer et al.l 
(200^. These two reactions (see, e.g.. IValle et al.ll2009l 
and references therein) affect both the extension and the 
morphology of the loop performed by intermediate-mass 
He-burning stars in the HR diagram. Hence the most 
recent values should always be adopted in stellar evo- 
lutionary codes in order to provide reliable theoretical 
predictions. 

Current calculations did not take into account He 
and metal diffusion, since its effect on the evolution- 
ary properties of intermediate-mass stars is minimal. 
To model super-adiabati c convection, we ad opted the 
mixing- length formalism (|Bohm-Vitensel[T958l ) . in which 
the efficiency of the convective transport depends on the 
mixing- length, i.e. I = aHp, where Hp is the pressure 
height-scale and a is a free parameter. The a value - 
1.74- was calibrated using the Standard Solar Models 
computed with both the same code and the same in- 
put physics, in particular using the same heavy-element 
mixture (As09), To account for current theoretical and 
empirical uncertainties, additional models with q;=1.9 
were also computed. We also account for convective core 
overshooting, during central hydrogen burning phases, 
by extending the central mixed region of lov—PHp from 
the border of the standard convective core defined by 
the Schwarzschild criterion, where /3 is a f ree parame- 
ter (|Castellani et al. II2000I : IValle et al.l l2009[ ). Note that 
we do not account for a mechanical core overshooting 
during the central He-burning phase. The use of the 
Schwarzschild criterion to determine the convective core 
extension, during these phases, would lead to an un- 
physical discontinuity in the radiative gradient at the 
border of the convective core. This effect is caused by 
the increase in the opacity due to the conversion of He 
into C and O (jCassisi et al.l 120041 ) . To overcome the 
problem we model the growth of the convective core 
an d the development of a s emiconvective region follow- 
ing [CaitellinriE^ (|1971al lH). Note that we neglect the 
convective overshooting at the base of the convective en- 
velope during shell hydrogen burning phases. 

2.2. Evolutionary tracks 

The approach we plan to adopt to estimate the in- 
trinsic parameters of the binary components does re- 
quire detailed set of evolutionary tracks and of isochrones 
properly covering the helium burning phases. The stel- 
lar masses of the evolutionary tracks range from 3 to 5 
M0, and the step in stellar mass is 0.05 Mq. There- 
fore, each set includes 42 evolutionary tracks. Un- 
fortunately, we still lack an accurate measurement of 
the iron and heavy element abundances of CEP0227. 
Therefore, we adopted the mean metallicity of LMC 
Cepheids ([Fe/H]=-0.33±0.13 dex) based on recent high- 
resolution, high sign al-to-noise spectra coll ected with 
UVES at ESO/VLT (|Romaniello et al.ll2008l) . The top, 
left panel of Fig. 1 shows the set of evolutionary mod- 
els computed at fixed helium (Y=0.260) and metal 
(Z=0.006) content and stellar masses ranging from 3 to 
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Fig. 1. — Top - Left - Hertzsprung- Russell Diagram of a set of evolutionary tracks computed at fixed metal (Z=0.006) and helium 
(Y=0.260) abundance and stellar masses ranging from 3 to 5 Mq. The plot is focussed around the evolutionary phases of central helium 
burning (blue loops). The free parameters adopted to deal vfith eonvective transport (mixing length, a), with conveetive core overshooting 
(/3) and with mass loss (rj) are also labeled. Top — Right - Same as the Top-Left, but for evolutionary tracks computed at fixed stellar mass 
(M=4.15 Mq), but different chemical compositions (see labeled values). Middle - Left — Same as the Top— Left, but for evolutionary tracks 
computed at fixed stellar mass (M=4.15 Mq) and chemical composition (Z=0.006, Y=0.260), but for two different assumptions concerning 
the mixing length (see labeled values). Middle — Right — Same as the Top-Left, but for evolutionary tracks computed at fixed stellar mass 
(M=4.15 Mq) and chemical composition, but for two different assumptions concerning the conveetive core overshooting. Bottom - Left - 
Same as the Top-Left, but for evolutionary tracks computed at fixed stellar mass (M=4.15 Mq) and chemical composition, but for three 
different assumptions concerning the mass loss rate. Bottom - Right - Same as the Bottom-Left, but for evolutionary tracks computed 
neglecting the conveetive core overshooting. 



5 Mq. The evolutionary phases plotted in this panel 
cover the central helium burning evolutionary phases 
(blue loops). 

To take into account possible changes in iron 
abund ance, due to intr in sic metallicity d isper- 
sion (IPompeia et all 120081: IRomaniello et all I2Q08I: 
IHunter et al.l l2009allH ) we also computed four 
independent sets of evolutionary tracks with 
metaUicity ranging from Z=0.004 to Z=0.008, 
while the helium content was fixed assuming 
an helium-to-metal en r ichment ratio AY / AZ=2 
(IPaeel fc Portinaril [l99l I J imenez et al.' '20031; iFlvnn! 
120041: IG ennaro. Prada Moron i. & Dcgl'lnnoccnti 2010) 
and a primo rdial helium abund an ce of 
(see, e.g., ICvburt et al.l [20Q1 ISteig mani 



y„= 0.2485 
I2Q06I: 



Peinibert. Luridiana. fc PeimbertI 120071 : iKomatiSu et alj 
2011[) . Therefore, current helium values range from 
Y=0.256 to Y=0.265 (see top, right panel of Fig. 1). 
We already mentioned that intermediate-mass stars are 



affected by uncertainties on the actual size of the He-core 
mass. To constrain the impact of the conveetive core 
overshooting during central hydrogen burning phases on 
the intrinsic parameters, the evolutionary tracks were 
computed either neglecting (/3o«=0) or accounting for a 
moderate conveetive core overshooting (/3o„=0.2, see top 
panels of Fig. 1). The two components of the binary 
system are relatively cool objects (see Table 1) with ex- 
tended conveetive envelopes. To constrain the depen- 
dence on the efficiency of conveetive transport in the 
super-adiabatic regions the calculations were performed 
assuming two different mixing length parameters. The 
former one -a=1.74- calibrate d using the standard solar 
model ()DeirOmodarme et al. 2011) and the latter one - 
a=1.90- to mimic a more efficient convection (see middle 
panels of Fig. 1). 

Recent theoretical investigations based on evolution- 
ary models accounting for both conveetive core over- 
shooting and for a pulsation driven mass loss suggested 
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Fig. 2. — Hertzsprung-Russell Diagram of sets of isochrones computed at fixed metal (Z=0.006) and helium (Y=0.260) abundance and 
ages ranging from 80 to 180 Myr, with a step in age of 10 Myr. 



that the latter physical mechanism can play a crucial 
role in settling the Cep heid mass discrepancy problem 
(jNeilson. Cantiello. fc L angcr 20"Tl[). A detailed analy- 
sis of the dependence of the blue loops on this physical 
mechanism is out of the aim of this investigation. How- 
ever, to constrain the dependence of the most probable 
solution on the mass loss rate, different sets of evolution- 
ary models were computed accounting for mass loss after 
the central hydrogen exhaustion. We used the Reimers 
formula 

dM/dt = x4 X 10-13 x L/gR [Mg/yr] 
where 77 is a free parameter, L,R and g are the lumi- 
nosity, the radius, and the surface gravity in solar units. 
Sets of models with four different 77 values, namely 0, 0.4, 
0.8 and 4.0, were computed. The evolutionary tracks 
constructed by assuming 7/=0.4 and 0.8 show marginal 
changes and in the following we will only take into ac- 
count the former ones (see bottom panels of Fig. 1). 

We ended up with a sample of more than 2,000 evolu- 
tionary tracks computed for this specific project. 



2.3. Stellar isochrones 



For each set of evolutionary tracks discussed in §2.2, we 
also computed a fine grid of stellar isochrones with ages 
ranging from a few tens to a few hundreds of Myr. To 
improve the accuracy of the parameter estimates (see 33]) 
and to provide robust determinations of the confidence 
interval, the isochrones were further interpolated. The 
actual grid of isochrones has a step in age of 0.02 Myr, 
necessary to cover the fast evolution of the stars across 
the blue loop. 

Subsets of isochrones for selected input parameters (see 
labeled values) in the HR diagram are shown in Fig l2.1l 

3. PARAMETER ESTIMATES 

The parameters that determine the evolutionary sta- 
tus of a star are its mass, /i, its age, r and its chemical 
composition, C,. The value of these parameters determine 
the position of a star in the HR diagram or, similarly, in 
any diagram were observable properties are displayed. 
Stellar models predict observables as a function of these 
three parameters, hence their values can be estimated 
by comparing predictions and observations. To provide 
robust estimates of both stellar masses and ages of the 
binary components we followed the Bayesian approach 
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described in lGennaro. Prada Moroni, fc Tognellil ()2012D . 
The interested reader is referred to the quoted paper, 
here we only summarize the relevant points of the ap- 
proach we adopted. 

3.1. Bayesian approach 

Given a set of observables q the probability of the 
model parameters, p = {ii, r, (^) is given by: 

cx C{p\q)f{p}- (1) 

Here C{p\q) is the Likelihood of the parameters p given 
the set of observations q and f{p) is the prior distribution 
of the parameters. 

We used three different pairs of observables, q, to 
determine the (/i, r, C) set of parameters of the binary 
components. In particular, we performed the compari- 
son in the q = (log L/i©, logTcff ), (M/Mq, i?/i?Q) and 
(logreff,log5) planes. 

To estimate the best values of the stellar masses and 
ages the probability f{p\q) is marginalized w.r.t the other 
two variables: 

G(r) = j £{t, C\q)f{T, fi, C) d^l dC (2) 
Hifi) - J £{t, f,, C|g)/(T, fi, C) dr dC (3) 

The mode of the marginal distributions is used as the 
best parameter estimato r (Jorgcnscn &: Linde gren 2005;; 
iGennaro. Prada Moron i, fc Tognellil l2012t) . ^Confidence 
intervals are defined by excluding 16% of the area under 
the posterior distribution on each side of the variable 
domains. 

In our treatment the metallicity is not considered as 
an unknown parameter to be determined. On the con- 
trary, we adopted the available prior information on 
the metallicity of Cepheids in the LMC to constrain 
our predictions on the stellar masses and ages. In 
detail, we adopted a Gaussian distribution of [Fe/H] 
with mean < [Fe/H] >= —0.33 and standard deviation 
"■[Fo/H] — 0-13 dex. The quoted values are based on the 
metallicity distrib ution of LMC Cepheids measured by 
iRomaniello et all (2008 ). 

Given the available information on the dynamical stel- 
lar masses we adopted two different priors for the stellar 
mass. In one case a flat prior is adopted, defined across 
the entire range of stellar masses for which we computed 
evolutionary tracks: M 6 [3,5]Mq. In the second case 
we used a prior with a Gaussian distribution with mean 
and standar d deviations based on th e dynamical mass 
estimates bv lPietrzvnski et all (|2010( ). 

Stellar ages were estimated separately for the two bi- 
nary components and for the system as a whole. The age 
of each star is determined from the G(t) marginal dis- 
tribution of equation ([2]). In the reasonable hypothesis 
that the two components of the binary system have the 
same age, their composite age distribution, i.e. the age 
of the system, can be simply derived as: 

Gc = Gp X Gs ; (4) 

where the subscripts indicate the composite (C), the 
primary component (P) and the secondary compo- 
nent (S) age distributions. In a recent investigation. 



IGennaro. Prada Moroni, fc Tognelhl ()2012D showed that 
the age distribution of a stellar system based on this ap- 
proach is a more precise indicator of the real age when 
compared to single stellar ages. 

The above approach was applied to each set of evolu- 
tionary models listed in Table 1, i.e. for each different 
assumption concerning the value of the overshooting pa- 
rameter -/3- the mass-loss -77- and the mixing length 
-a- parameter. The posterior probabilities were com- 
puted using both the flat and the Gaussian prior on the 
stellar mass. Selected results of the individual solutions 
are given in Table 3. The solutions based on different 
observables are discussed in the following subsections. 

3.2. The Hertzsprung-Russell Diagram 

As a first step to constrain the intrinsic parameters 
of the binary system, we applied our method in the 
Hertzsprung-Russell Diagram. The main advantage in 
using this plane is that we can estimate the intrinsic lumi- 
nosities using the Stefan-Boltzmann relation. Therefore, 
the solution is independent of distance modulus, redden- 
ing and stellar mass. The main drawback is that the two 
parameters (luminosity and temperature) are correlated 
and more affected by measurement errors than the stellar 
mass and the radius (see Fig. 3). However, our method 
takes into account the fact that temperature and lumi- 
nosity are not independent, and indeed the covariance 
matrix is also included in the solution (see equations (3) , 
(4) and (5) in lGennaro. Prada Moroni, fc Tognellill2012| i. 

3.3. The stellar mass vs radius plane 

To further constrain the intrinsic parameters of the 
binary system we also applied the Bayesian method in 
the stellar mass vs radius plane. The main advantage 
in using this plane is that the adopted observables are 
independent of distance modulus and reddening and they 
are also independent from each other (see Fig. 4). 

3.4. The logTeff vs logg plane 

Finally, we applied the Bayesian method also into the 
effective temperature vs surface gravity plane. The main 
advantage in using this plane is that the adopted ob- 
servables are independent of distance modulus and red- 
dening. Moreover, the observables are independent from 
each other and the isochrones display a relevant differ- 
ence in this age range (see Fig. 5). This plane com- 
bines three independent measurements, thus providing 
the most stringent test for the theory. 

4. CONSTRAINTS ON THE STELLAR MASS AND 
AGE OF THE BINARY SYSTEM 

The most probable mass and age values that we found 
indicate that the mixing length parameter has a minimal 
impact on the estimate of these parameters. This result 
is not surprising, since the extent in effective tempera- 
ture of the blue loop is minimally affected by the change 
in a from 1.74 to 1.90. This is the reason why in Ta- 
ble 3 we only included the results based on models with 
a=1.74. What is clear from the values listed in Table 3 
is that the solutions based on evolutionary models that 
do not account for the convective core overshooting and 
those with enhanced mass-loss (77 = 4.0) are in worse 
agreement with the measured dynamical masses, when 
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Fig. 3. — Top - Left- Comparison in the Hertzsprung-Russell Diagram between the binary components and the isochrone giving the 
maximum of Gq = Gp X Gg (see Table 3). In this panel the solution is based on evolutionary models constructed by neglecting convective 
core overshooting (/3=0) and mass loss (?7=0). The solid and the dashed lines show the isochrones with most probable age derived using 
the Gaussian and the flat prior on the mass distribution, respectively. The blue bar marks the position of the Cepheid (primary), while the 
red ones the position of the red giant (secondary). The intrinsic parameters of the isochrones are also labeled. Top -Right- Same as top 
left, but for models accounting for mild (/3=0.2) convective core overshooting. Bottom - Same as the top right, but for models accounting 
for a canonical (left, r)=QA) and an enhanced (right, »7=4.0) mass loss rate. 



compared with the results based on models accounting 
for the convective core overshooting and 77 = 0.0 or 0.4. 

In the case of a flat mass prior, the /3 = 0.0 models 
predict stellar masses that are systematically larger than 
observed, even though the difference with the dynamical 
mass is typically within the la confidence interval. Sim- 
ilarly, the models with enhanced mass-loss, in the case 
of flat mass prior, tend to underestimate the mass of the 
primary component. These two sets of models give a 
good agreement with the data, only when the method is 
used in the M-R plane. 

Comparing the results for a flat and a Gaussian mass 
prior, it is clear that the latter drives the solution to- 
wards the measured dynamical mass. While this might 
appear obvious, we will show in Sect. 14.11 that not only 
the mode of the mass distributions are shifted towards 
the measured values, but also the whole solution has a 
higher reliability, measured using the so-called evidence 
parameter (see §4.1). The impact of the Gaussian prior is 



very important also for the age results. In the Gaussian 
mass prior case, the age solutions for the two components 
are generally more precise, meaning that the uncertainty 
interval is narrower compared to the flat prior case. In 
addition to this, there is also a better agreement between 
the ages of the two components, and therefore, a more 
robust result for the composite age of the system. 

For each set of models, we plotted the isochrones corre- 
sponding to the most probable composite age (indicated 
by a C in Table 3) in Figs. 3, 4 and 5. The three figures 
show the most probable isochrone as obtained by apply- 
ing the method in the L-T, in the M-R and in the g-T 
plane. Solid and dashed lines display the most probable 
ages for a flat and a Gaussian mass prior, respectively. 

In passing, we note that evolutionary models con- 
structed assuming a slightly enhanced mass loss rate 
(77=0.8) also provide very similar solutions to the 
canonical (77=0.4) case. This indicates that canoni- 
cal recipes to account for mass loss rate do not al- 
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Fig. 4. — Top —Left- Comparison in the stellar mass vs radius plane between the binary components and the isochrone giving the 
maximum of Gc = Gp X Gg (see Table 3). The symbols and the lines are the same as in Fig. 3. Top -Right- Same as top left, but for 
models accounting for mild (13=0.2) convective core overshooting. Note that the dashed line overlaps with the solid line, since they attain 
very similar values in the M-R plane. Bottom - Same as the top right, but for models accounting for a canonical {r]=OA, left) and an 
enhanced (»7=4.0, right) mass loss rate. 



low US to provide more quantitative constraints on 
mass loss efficiency during advanced evolutionary phases 
(jNeilson. Cantiello. fc Langeii 120111) . However, a de- 
tailed treatment of the pulsation-driven mass loss inside 
the instability strip is out of the aim of this investigation. 

4.1. Evidence of the most probable solutions 

The confidence intervals of the different solutions pro- 
vide a hint concerning the robustness of individual fits 
using the different observables. However the precision 
of the result alone migth be a deceptive indicator of the 
goodness of the fit between theory and data. A nar- 
rower confidence interval in the (/x, r) parameter space 
does not necessarily imply that the corresponding choice 
either of the observables, or of the mass prior or of the 
set of models is better than other possible choices. 

A quantitative measurement of the relative goodness of 
the different solutions can be provided by the Bayes Fac- 
tor. The Bayes Factor, BFy, between two sets of models 
is defined as the ratio of the evidence of the two mod- 
els. The evidence itself is the integral of the Likelihood, 



marginalized o ver the model para meters prior distribu- 
tion (see, e.g., lBailer-Jonesll20lj]) . The evidence gives 
the total probability of the observed data given the par- 
ticular set of models considered. Since our sets of models 
differ for the choice of the mixing length parameter, a, 
the core overshooting parameter, /3, and the mass-loss 
parameter, rj, we introduce = {ai,f3i,r]i) to identify a 
particular choice of the three parameters. The evidence 
is therefore: 

£M EE 8{q]Ei) = J /:(r,Ai,C|<r,S)/(T,M,C|S)drdAidC 

The evidence depends on the plane in which the com- 
parison between data and models is performed (L-T, 
M-R, g-T), on the choice of the prior distribution of 
the model parameters (i.e., /(r, ^|S), in our case flat 
or Gaussian mass distribution and Gaussian metallicity 
distribution) and on the S-set of parameters identifying 
that set of models (see Fig. 6). £i{q) is proportional to 
the total probability of observing the data q under the 
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Fig. 5. — Top -Left- Comparison in the logarithmic effective temperature vs logarithmic surface gravity plane between the binary 
components and the isochrone giving the maximum of Gq = Gp X Gg (see Table 3). The symbols and the lines are the same as in Fig. 3. 
Top -Right- Same as top left, but for models accounting for mild (/3=0.2) convective core overshooting. Bottom - Same as the top right, 
but for models accounting for a canonical (left, »y=0.4) and an enhanced (right, r)=4.0) mass loss rate. 



assumption that the 5^ model is the correct model. 

In equation ([5]) the q observables refer to a single star. 
The composite evidence for both stars is simply: 

£i{{qp, qs}) = £t{qp) x £i{qs) (6) 

From equation ([5]), it follows that the Bayes factor is: 



BF, 



(7) 



The Bayes factor does not constrain the triplet (r, /x, Q 
that provides the highest posterior probability for each 
set of models. However, it can be used to quantify which 
set of models, variable plane and prior distribution choice 
give the best overall agreement with the data. Typically, 
in order to claim that a set is preferred over another, 
the Bayes factor should be s maller /larger than one , with 
either BF < 0.1 or BF > 10 (jKass fc RftervlfTOOl . 

Table 4 gives the BF values using as a reference the set 
of models with a — 1.74, j3 — 0.2 and ?/ = 0.4, together 
with the variables q — (log 5, logTcff) and with the Gaus- 
sian mass prior distribution. This reference set gives the 



largest evidence for the composite system, therefore it 
can be considered -within our range of models- the best 
set of models and the best prior distribution choice to de- 
scribe the observations. The BFs listed in Table 4 were 
estimated by dividing the corresponding evidence by the 
evidence of the reference set either for the primary star 
(P entries), or for the secondary star (S entries) or for 
the composite evidence (C entries). 

The (/3 = 0.2, -q = 0.0) set of models gives BFs for the 
primary (see Fig. 6), the secondary (see Fig. 7) and the 
composite system that are very similar to the reference 
set. Therefore, these two sets give similar likely solutions 
and with the present data we can not distinguish between 
them. On the contrary, the set without overshooting, 
/3 = 0, and the set with enhanced mass-loss rate, 77 — 
4.0, show very low values of the BFs, suggesting that 
these parameters can be excluded, hence, confirming the 
qualitative analysis at the beginning of this section. 

It is worth noticing that, the BFs for a given q plane 
and a given S set are generally larger when the Gaussian 
mass prior is used. This suggests that the likelihoods 
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themselves are already centered in a region of the pa- 
rameter space which is close to the dynamical mass mea- 
surements. Therefore, the Gaussian prior enhances the 
posterior probabilities when compared with the flat case. 
This finding is also supported by the fact that several so- 
lutions obtained using a flat mass prior distribution are 
characterized by large confidence intervals (see top pan- 
els in Figs. 6 and 7 and the values listed in Table 3). In 
addition to this, the plane in which the largest BFs are 
found is the g-T plane. This means that the models, 
when the measurements of stellar mass, radius and ef- 
fective temperature are combined, are still able to repro- 
duce these three independent observables, and actually 
the total probability increases. 

Finally, we mention that the presence of two different 
regions in all the probability diagrams for the secondary 
component (see Fig. 7) is the consequence of an intrinsic 
evolutionary feature during the blue loop (helium burn- 
ing phases). Intermediate- mass stars soon after helium 
ignition (tip of the red giant branch) move, at fixed stel- 
lar mass, towards higher effective temperatures (the so- 
called blue tip) and then back to lower effective temper- 
atures (jBono et al.ll200dl and references therein). Dur- 
ing these phases stellar structures with similar masses 
can be found in similar positions at different ages. How- 
ever, the contribution of the two probability-regions to 
the marginal distribution is very different, with the main 
region contributing a much larger probability than the 
secondary one. The latter region only causes an increase 
in the skewness of the marginal distributions, without 
manifestation of secondary maxima. 

5. THE COLOR-MAGNITUDE DIAGRAM 

The binary system CEP0227 offers a unique oppor- 
tunity to test the plausibility of the physical assump- 
tions adopted in evolutionary and pulsation models be- 
cause the comparison between theory and observations 
can be performed by using observables (luminosity, ra- 
dius, effective temperature, surface gravity) directly pre- 
dicted by evolutionary models. This means that these 
observables do not need to be transformed into the ob- 
servational plane. However, the Cepheid mass discrep- 
ancy problem arises from the comparison between evolu- 
tionary prescriptions and observations in the CMD. The 
main advantage in using this plane to estimate the evo- 
lutionary mass is that the estimate only relies on two 
mean magnitudes in different photometric bands. The 
main drawback is that the mass estimates are affected 
by uncertainties in the true distance modulus and in the 
reddening. To overcome the latter problem the mass es- 
timates have also been performed in CMDs based on NIR 
bands (Bono ct al. 2001). 

Absolute distance determination and NIR magnitudes 
for the binary CEP0227 are not available yet. We have 
already collected new NIR time series data and we plan 
to provide a detailed analysis of the CMD and of the 
color-color plane in a forthcoming investigations. 

However, the solutions based on theoretical observables 
do allow us to constrain the possible occurrence of thorny 
systematic uncertainties in the bolometric corrections 
and in the color-temperature relations adopted to trans- 
form theoretical predictions into the observational plane. 
Moreover, we can also estimate the reddening of the sys- 
tem by assuming as a prior a Gaussian distribution on the 



true distance modulus. We used the theoretical solution 
which gives the largest evidence for the composite sys- 
tem, i.e. the 151 Myr isochrone for the = 0.2, rj = 0.4 
set. This is the most probable age for the composite 
G{t) in the case of the g-T plane. We transformed the 
isochrone into the Johnson-Kron-C ousins photometric 
system using the spectra provided bv lCastelli fc Kuruczl 
()2003() based on the ATLAS9 model atmospheres (see 
Fig. 8). We also performed a test by using the spec- 
tra based on the atmosphere mo dels pr ovided by the 
PHOENIX code in the version by iBrott fc Hauschildt] 
(200i>) ■ We found that the difference -in the typical tem- 
perature, metallicity and surface gravity regime of the 
binary we are dealing with- is negligible, therefore in the 
following we adopt the evolutionary models transformed 
using the ATLAS9 models. 

The comparison into the observational plane (CMD) 
between the isochrone of the most probable solution and 
observations provides independent estimates of both the 
reddening and the true distance modulus of the binary 
system. The reddening in t he V ,I-bands was fixed by 
adopting the iCardelli et al.l (|1989[ ) extinction law. The 
solution was found by maximizing the probability: 



p{E{B -V),DM)^ Yl 



1 



^=PS 27rcry,cr(v_/)^ 



X (8) 



exp 



Xv^ + X(v-i), {DM - DMoY 



dm 



where: 



XV,^ 



V^ - V{m; E{B - V), DM) 



(TV, 



and 



(y ~ I), -{V- I){m- E{B - V), DM) 



The index i runs on the primary and on the secondary 
star (P, S), Vi and {V — I)i are the observed mean mag- 
nitudes and colors with their errors cry. and (T(y_/) . . The 
theoretical quantities are function of the mass along the 
isochrone and of the adopted distance modulus and red- 
dening. The integral is performed along the isochrone 
integrating over the mass of the models. Given the avail- 
able information on the LMC distance modulus, we used 
a Gaussian prior with mean and standard devi ations 
given by /X = 18.45 ± 0.04 mag (i Storm et ani201lD . 

Fig. 8 shows the 151 Myr isochrone in the V, V-I 
CMD, with a reddening of 0.142t°;o"g mag and a dis- 
tance of 18.531q'q2 niag as foimd by maximizing equa- 
tion ([U . This estimate of the reddening agrees very well 
with the reddening estim ates based on Schlegel's map, 
name ly 0.15±0.03 mag (jSchlegel. Finkbeiner. &: Davis! 
|1998[ ). The same outcome applies to the LMC distance, 
since it agrees with similar estimates available in the lit- 
erature (jStorm et al.ll20lTl and references therein). This 
evidence is also supported by the fact that the plane of 
the LMC is also tilted along the line of sight and we 
are not applyin g any correction for the position of the 
binary system ijPersson et al.l l2004f) . In order to con- 
strain the sensitivity of both distance modulus and red- 
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Fig. 6. — Posterior probability distributions for the primary component according to the set of evolutionary models constructed by 
assuming /3 = 0.2, rj = 0.4 in the L-T (left), in the M-R (middle) and in the g-T (right) plane using either a flat (top panels) or a 
Gaussian (bottom panels) prior on the mass distribution. Each 2D probability is obtained after marginalizing over the metallicity Gaussian 
distribution for LMC Cepheids. In each panel the probability was normalized to its maximum value for displaying purposes, therefore the 
color coding range from (blue) to 1 (red) normalized probability. The top and the right insets display the marginalization for the stellar 
age and the stellar mass of the primary, respectively. The vertical dashed lines show the confidence intervals. The y-axis of the insets show 
the actual value of the marginal distributions before the normalization. 




Fig. 7. — Same as Fig. 6, but for the secondary component. 
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Fig. 8. — The isochrone of the most probable solution (13=0.2, 
r]=OA, g—T plane, Gaussian mass prior) transformed into the V,V- 
I color-magnitude diagram. The blue and the red symbols display 
the Primary (Cepheid) and the Secondary component of the binary 
system. The true distance modulus and the reddening estimated 
using the combined system are also labeled. 

dening on the prior we also adopted a flat prior between 
^=18 and 19 mag. Note that the adopted interval is 
significantly larger than the range of LMC distances in- 
dicated in literature and we found a true distance mod- 
ulus of ^—IS.blto'^l mag and a reddening of E(B-V) = 
0.12751q'o7^^ mag. The new estimates have larger confi- 
dence intervals, but they agree quite well with the esti- 
mates based on the Gaussian prior. The above findings 
and the robustness of the most probable solution (sec 
Fig. 9) provide an independent support to the current 
set of atmosphere models a dopted to transform theo ry 
into the observational plane (|Castelli fc Kuruc^l2003[ ). 

6. SUMMARY AND FINAL REMARKS 

The discovery by OGLEIII of eclipsing binary systems 
in the Large Magellanic Cloud hosting classical Cepheids 
provided the unique opportunity of measure their masses 
with a precision ranging from 1% (CEP0227) to 1.5% 
(CEP1812, paper II). It goes without saying that the new 
measurements will allow us to constrain the precision 
and the plausibility of the the input physics and physical 
assumptions adopted in constructing both evolutionary 
and pulsation models. In this investigation, we addressed 
the evolutionary mass and the age of the binary system 
CEP0227. We adapted the general Bayesian approach 
developed bv lGennaro. Prada Moroni, fc Tognellil (|2012[ ) 
to constrain the physical parameters of the binary com- 
ponents. 

In order to provide robust statistical solutions, we com- 
puted several sets of evolutionary models by covering a 
broad range of chemical compositions (see Table 1) and 
stellar masses. To constrain the dependence on the free 
parameters currently adopted in canonical evolutionary 
models, independent sets of models were constructed by 
assuming two different values of the mixing length pa- 
rameter (q;=1.74, 1.90). Moreover, to account for the 
possible occurrence of extra-mixing during central hy- 
drogen burning phases, the models were constructed ei- 
ther by neglecting or by including a moderate convective 
core overshooting (/3ou=0.2). Finally, to account for the 
mass loss during evolved evolutionary phases the mod- 
els were constructed either by neglecting or by including 
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Fig. 9. — Robustness of the most probable solutions for the true 
distance modulus and the reddening of the binary system obtained 
using the isochrone with the largest evidence and a prior on the 
true distance modulus (18.45±0.04) mag. The color coding shows 
the transition from less (blue) to more (red) accurate most probable 
solutions. 

mass loss. To mimic different efficiencies in the mass loss 
rate, we constructed sets of models accounting either for 
canonical (7^=0.4,0.8) or for enhanced {rj—A) mass loss 
rates (see Table 1). 

To overcome possible deceptive errors in the estimate 
of the physical parameters, the comparison between 
theory and observations was performed in three differ- 
ent planes: luminosity vs effective temperature (logL, 
logTe//), mass vs radius (M,R) and surface gravity vs 
effective temperature (logg, logTe//). We computed an 
independent solution (mass, age, chemical composition) 
in the three quoted planes for every set of evolutionary 
models we have constructed. Note that individual so- 
lutions were computed by adopting two different priors 
for the stellar mass. We adopted either a fiat distri- 
bution over the entire range of stellar masses covered 
by computations or a Gaussian distribution anchored to 
the mean and the standard deviation of the dynamical 
mass measurements (paper I). We also adopted a Gaus- 
sian prior for the metallicity, in particular we adopted 
the mean and the sta ndard deviation fo r LMC Cepheids 
recently provided by iRomaniello et al.l ()2008D by using 
high-resolution, high signal to noise spectra. 

The main results we obtained in the comparison be- 
tween theory and observations are the following; 

• The most probable solutions are minimally affected 
by the adopted mixing length parameter, within 
the range covered by our models. 

• The most probable solutions based on a Gaussian 
mass prior are typically more accurate than those 
based on the flat mass prior. 
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• The most stringent test to constrain the models 
were obtained in the g-T plane, due to the fact 
that the adopted parameters involve three indepen- 
dent observables, namely mass, radius and temper- 
ature. In this plane, the most probable solutions 
were obtained by using the evolutionary models 
that account for a moderate convective overshoot- 
ing (/3o„=0.2) and canonical mass loss rate (77=0.4). 
In particular, we found that the mass of the pri- 
mary (Cepheid) is M=4.14l°;°^ Mg, while for the 
secondary is M=4.15lo;o5 ^0- ^o^e that current 
estimates do agree at the 1% level with dynam- 
ical measurements. We found ages for the indi- 
vidual components and for the combined system 
(t=15l3 Myr) do agree at the 5% level, support- 
ing once again the robustness of the solution. It is 
noteworthy that similar solutions are obtained by 
using moderate convective overshooting (/3o„=0.2) 
and no mass loss (ry=0). This indicates that ei- 
ther a vanishing or a canonical mass loss rate 
play a marginal role in the Cepheid mass discrep- 
ancy. Current findings support the preliminary re- 
sults concerning the evolutionary mass of CEP0227 
provided by Cassisi & Salaris (2011). They also 
suggest that evolutionary masses based on evolu- 
tionary models accounting for moderate convec- 
tive core overshooting -a typical assumption in 
the comparison between theory and observations 
([Pietrinferni et al.|[2006l) - agree quite well with the 
dynamical mass of Cep0227. The overall agree- 
ment among dynamical, evolutionary and pulsation 
masses will be addressed in a forthcoming paper in 
which we plan to constrain the pulsation mass us- 
ing independent observables. 

• The most probable solutions based on the mod- 
els constructed by assuming a moderate convec- 
tive overshooting (/3ot)=0.2) and an enhanced mass 
loss rate (^7=4) present large confidence inter- 
val when compared with the solutions based on 
canonical mass loss rate (see Table 3). Note 
that current finding should be cautiously treated. 
The reason is twofold. The difference be- 
tween the estimated masses and the dynamical 
masses, taken at face value, do agree at the 1% 
level. We are assuming a very crude approxi- 
mation to account for the mass loss, while cur- 
rent empirical and theoretical findings indicate 
that the mass l oss might be driven by the pulsa- 



tion instability 


Neilson. Cantiello. & LanKeiil2011l: 


IMatthews et al. 


120111). This hypothesis is fur- 



thcr supporte d by t he recent results obtained by 
iNcilson et al. (|2011[ ). They found that the ob- 
served period change of Polaris implies a mass 
loss rate of 10~^ Mq yr~^, and in turn a new 
firm empirical argument for enhanced mass loss in 
Cepheids. 

The above discussion concerning the formal precision 
of the results relies on the width of the confidence in- 
terval. However, a solution characterized by a narrow 
confidence interval might not be the most probable so- 
lution. The Bayes Factor between two sets of models 



is a more robust statistical parameter to constrain on a 
quantitative basis the most probable solutions. 

By using the Bayes Factor the main results we obtained 
in the comparison between theory and observations are 
the following: 

• The solutions with the largest evidence were ob- 
tained in the g-T plane by using sets of models 
that account for a moderate convective core over- 
shooting (/3ou=0.2) and either a vanishing (77=0) or 
a canonical mass loss rate (77=0. 4) give very similar 
results. This finding supports the results based on 
the confidence interval. 

• The solutions based on sets of models that cither 
neglect convective core overshooting or account for 
an enhanced mass loss rate (ti—A) are characterized 
by very low Bayes Factors. Thus suggesting that 
these two sets of models give a poorer description 
of the observed properties of the binary system. 

We also transformed the isochrone of the most 
probable solution {(3=0.2, 77=0.4, g-T plane, Gaus- 
sian mass prior) mentioned above into the observa- 
tional plane using the at mosphere models provided by 
ICastelh fc KTinica ((20031 ). By assuming a Gaussian 
metallicity prior distribution for the Cepheid metallicity 
and by assuming the Cardelli et al. (1989) reddening law, 
we found a true distance modulus -18.53lg'Q2 ™&g~ ^ind 
a reddening value -E(B-V)= O.U2toml mag- that agree 
quite well with similar estimates available in the litera- 
ture ([Schlcgel, Finkbcincr, & Davis 1998; Pcrsson et al] 
12004 Istorm et all 120 111 ). The use of a fiat prior on 
the LMC distance, minimally affects these results. A 
more detailed discussion concerning new independent 
distances of the binary system and of the Cepheid will 
be addressed in a forthcoming investigation. 

Intrinsic variables in eclipsing binary systems are a fun- 
damental laboratory to constrain the input physics of 
both evolutionary and pulsation models. This was con- 
sidered till a few years ago a broken ground for quanti- 
tative astrophysics, due to the limited number of known 
systems and the lack of homogeneous and precise pho- 
tometric and spectroscopic data. The unprecedented 
photometric catalogs released by microlensing experi- 
ments (EROS, OGLEIII) tog ether wit h currently under- 
way optical fPanSTA RRS-l.lStubbs e t al. (2010); Palo- 
mar Transient Factory , Ivan Evken et al.i (|2011|): CHASE 
iFoerster. et all (|20Tol )) and NIR rVVV. iMinniti et alj 
(|2010D ) surveys and the spectrographs available at 8m 
class telescopes, are a very good viaticum in waiting for 
GAIA. 
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TABLE 1 

Intrinsic stellar parameters adopted to compute evolutionary tracks. 
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^ Metal (Z) and helium (Y) abundance by mass. 
^ Mixing length parameter. 
^ Overshooting parameter. 
Mass loss parameter. 



TABLE 2 

AnorTEn iktriksic rAnAiiETKRS for the iiixar^- s^'stei[ rEP0227. 



Parameter Value Ref.* 



Mass (M/Mq)^ 


4.14±0.05 






4.14±0.07 




Radius (R/Rq)^ 


32.4±1.5 






44.9±1.5 






5,900±250 




5,080±270 




V(mag)''''^ 


15.940±0.001 




16.06±0.02 




/(mag)'"''^ 


15.246±0.001 




15.11±0.01 




fj, (mag)'* 


18.45±0.04 


2 


E{B - V)(mag)° 


0.15±0.03 


3 


[Fc/ll] (dcyiY 


-0..3.3±0.05 


4 



^ References: 1) Pietrzynski et al. (2010); 2) Storm et al. (2011); 3) Schlegel et al. (1998); 4) Romaniello et al. (2008). 
^ The first line, for each parameter, refers to the primary (Cepheid), while the second one to the secondary (companion). 

Apparent V and I-band magnitudes. The mean magnitude of the Cepheid was estimated in flux and then transformed into magnitude. 

True distance modulus. 

® Reddening according to the Schlegel's map. 
f 

Mean iron abundance (scaled-solar mixture). 



Valle, G.; Marconi, M.; Degl'Innocenti, S.; Prada Moroni, P. G. 

2009, A&A, 507, 1541 
van Eyken, J. C. 2011, AJ, 142, 60 



Venn, K. A. 1999, ApJ, 518, 405 
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TABLE 3 

Most probable ages and masses for the different observables, the different sets of evolutionary models and for the 

PRIORS adopted for THE STELLAR MASS DISTRIBUTION. 

13 = 0.0, ri = /3 = 0.2,T) = /3 = 0.2, 77 = 0.4 /3 = 0.2, r? = 4.0 
MP"! Obs.'' Star'^ A.gc<* Mass"* Age'= Mass** Age*' Mass* Agc^ MassS 
[Myr] [Mrs,] [Myr] [Mq] [Myr] [Mq] [Myr] [Mq] 

P 115ir 4.50l«:- 1651?;^ 3.98l»i^ 167l^« 3.97^^^ 137+^, 3.601;^^ 

L-T s 96+_f 4.70l°i« 143l?« 4.15i°;i^ UT+H 4.15+',f, U8+_ll ^.W+'ofs 

C W7t\l IS^l^i 150+1'^ 137+^^ 

P 134+" 4 13+0-04 154+4 . 1Q+0.07 153+6 4 14+0.04 145+21 4 14+0.05 

M-R S 114tf 4.19t°;°| 1451^ 4.14l0:0| 14611 4.15t°;°| Udtl 4.15t°;°f 

C 144t3 150+5 I50t| 1451J5 



Flat 



P 115t« 4.50l«:^;; 16611? 3.981";;;^ 

g-T S 88lf ^-^St'oM 1431*} 4.15l^;2« 

C 95+" 154+2? 149+f 139+^« 

-3 -14 -9 _S2 



WYtll 3.96l«:?4 226lf^ 3.45lEJ:»° 

1y|P;+37 /I 1^+0.21 1/11+36 . nn+O-lS 



h 



p 144+5 4 -,0+0.05 -.r.+r, 4 ^o+O-OS ,54+4 4 ,,+0.04 ,07+6 4 15+0.04 

L-T S llSif 4.171°:°? 1471^ 4.151°;°^ 148t« 4.151°;°^ 150l^ 4.Ut°f, 

C 145l5 IS^ig ISllg 1411^1 

p 134+11 4,4+0.02 154+3 4,4+0.03 i„+4 4 14+0.03 145+17 4 14+O.03 

r IJ^-e ^•l''_o.04 10''_4 ^•1^_0.04 J^O'J-a ^-^^-0.04 ^_7h ^-^^-0.04 

Gauss '^-I^ S 148+^'' 4 17+0-04 ,45+6 4 ,4+0.04 ,4^+6 . 15+0.02 ,40+7 4 -,4+0.05 

O l^o_34 ^•li_o.07 l^O-S *-l^-0.04 ^™-3 ^-^^-0.05 l™-3 ^-1^-0.04 



C 1461^ 150l^ 



1+4 



P 145t| 4.13lH^ 154i! 4.1ll0;05 154+4 4.14+0.04 ,37+19 415+0.04 

g-T S 148l|2 4.17l°;°f 14711 4.15lH| 147t! 4.151°;°^ 1491^ 4.14; 



14611 



15313 151+4 



+0.07 
0.05 



3 



Prior on the stellar mass distribution, either flat or Gaussian distribution. 
^ Observables adopted in the solution to constrain the age and the mass of the binary system: L— T: luminosity and temperature; M-R: stellar mass and radius; g— T: 
surface gravity and temperature. 

^ Solution for the primary (P, Cepheid), secondary (S, red giant), combined (C). 

^ Solution for stellar age (Myr) and stellar mass (solar units) based on evolutionary models with an overshooting parameter ^=0 (no overshooting) and a mass— loss 
parameter 77 = (no mass-loss). 

Solution biisod on evolutionary models with p = 0.2 (mild overshooting) and j] = (no mass-loss). 
* Solution based on evolutionary models with /3 — 0.2 (mild overshooting) and 77 — 0.4 (canonical mass-loss rate). 
^ Solution based on evolutionary models with fB=0.2 (mild overshooting) and t;=4.0 (enhanced mass-loss rate). 

Errors in italic indicate that the confidence interval is poorly defined. 
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TABLE 4 

Bayes Factors for the different sets of models and for either a flat or a Gaussian mass prior. 



MP'' Obs> 


Star': 


= O.Q, ri = 0"^ 


13 = 0.2, r? = 0" 


/3 = 0.2, r] = 0.4f 


/3 = 0.2, )7 = 4.0S 




P 


0.395 


0.393 


0.369 


0.007 


L-T 


S 


0.367 


0.247 


0.255 


0.248 




c 


0.145 


0.097 


0.094 


0.002 




p 


0.014 
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< 10-3 


Flat M-R 
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< 10-3 


0.001 
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< 10-3 


< 10-3 


< 10-3 


< 10-3 
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0.011 
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0.311 
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0.003 



Gauss 
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1.000 


0.978 
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Prior on the stellar mass distribution, either flat or Gaussian distribution. 
^ Observables adopted in the solution to constrain the age and the mass of the binary system: L— T: luminosity and temperature; M— R: stellar mass and radius; g— T: 
surface gravity and temperature. 

Solution for the primary (P, Cepheid), secondary (S, red giant), combined (C). 
^ Solution for stellar age (Myr) and stellar mass (solar units) based on evolutionary models with an overshooting parameter ^=0 (no overshooting) and a mass— loss 
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